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0\ ', Abstract 

, The Generalized Lotka Voltera (GLV) formalism has been introduced in order to 

explain the power law distributions in the individual wealth {wi{t)) (Pareto law) 
and financial markets returns (fluctuations) (r) as a result of the auto-catalytic 
, (multiplicative random) character of the individual capital dynamics. 

^ ■ As long as the multiplicative random factor (A) is extracted from the same prob- 

ability distribution for all the individuals, the exponent of the power laws turns out 
I to be independent on the time variations of the average (< A >). This explains also 

O ■ the stability over the past century of experimentally measured Pareto exponent. 



In contrast to the scaling properties of the single time ("unconditional") probabil- 
ity distributions, the (auto-)correlations between observables measured at different 
times are not correctly reproduced by the original GLV, if the variance (cj^) of A is 
' time independent. In the GLV formalism the volatility (r^) auto-correlations decay 

exponentially while the measurements in real markets indicate a power law with a 
very small exponent. 

We show in the present paper that by making the variance of the individual wealth 
changes cr^ a function of the market volatility < >, one correctly reproduces the 
market volatility long range correlations. 

Moreover, we show that this non-trivial feedback loop between the market price 
volatility and the variance of the investors wealth leads to non-trivial patterns in 
the overall market trends. If the feedback is too strong, it may even endanger the 
market stability. 

Key words: volatility auto-correlations, power laws, econophysics, Lotka Volterra, 
stochastic logistic, behavioral finance 
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1 Introduction 



In the last decade, the Microscopic Representation techniques were used in a 
wide range of subjects[l]. In particular, Levy Levy and Solomon [2,3] (LLS) 
have devised a model of the financial markets in terms of a large number of 
virtual investors characterized each by a current wealth, portfolio structure, 
probability expectations and risk taking preferences (for a review see [4]). Such 
models allowed to uncover and study: 

• Market effects of arbitrarily inhomogenous and non-rational traders behav- 
ior. 

• Returns stochastic properties: autocorrelations, volatility, trading volume 

• Predation, competition and symbiosis among species. 

• Heavy-tailed market returns distributions related to the ratio between the 
capital entering the market and the increase in market stock capitalization. 

In order to understand what are the crucial factors governing this complex 
dynamics displayed by the microscopic representation of markets, one con- 
structed more schematic models which, while discarding some of the realistic 
features of LLS, still conserve the crucial dynamical features of the market. 

The GLV is such a model that embodies some styhzed features of the LLS 

model in a more generic framework. Instead of following in detail the way the 
market price influences each investor population and individual i , it was as- 
sumed that this influence can be represented through multiplying their wealth 
Wi{t) by stochastic multiplicative factors Xi{t). This is naturally suggested by 
the LLS model simulations in which the investments of the individuals (and 
consequently their returns) are fractions of their wealth (as implied by the 
constant relative risk aversion utility functions). This is also consistent with 
the recent measurements by [5] of the exponent of the power law distribution 
of the market order volumes. 

The stochastic proportionality between personal returns and personal wealth 
is consistent with the real data that show that the (annual) individual income 
distribution is proportional to the individual wealth distribution [6] . We pro- 
posed [7-11] therefore a model including the above stochastic autocatalytic 
properties of the capital as well as the cooperative, diffusive and competitive/ 
predatory interactions between the investors. The GLV described below is a 
straightforward stochastic generalization of the Lotka-Volterra system (and 
of the discrete logistic equation) well known previously in population biology 
and social sciences. 

As explained below, it automatically leads to many of the well known ex- 
perimental features of the real markets. However, some of the initial GLV 
simplifications were too drastic. In particular the assumption that the indi- 
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vidual returns are extracted from a probability distribution with fixed variance 
lead to the result that the market volatility auto-correlations decays exponen- 
tially with time. This is in stark contradiction with the measured real market 
properties [12-16]. In the present paper we identify the microscopic dynamical 
features which are responsible for "volatihty clustering" effect: the fact that 
the variance of the individual invested wealth changes is influenced by the 
global market price volatility. 

We show that this leads also to the emergence of a feedback loop which may 
in certain conditions destabilize the market. 



2 Background on the simple GLV model 

More than a hundred years ago, Pareto [17] discovered that the number of 
individuals with wealth (or incomes) with a certain value w is proportional to 

It turns out that in the conditions in which the participants in the market do 

not have a systematic advantage one over the other (which is in fact expected 
in an efficient market), realistic market dynamics of the LLS and GLV types 
lead always to Pareto laws. 

Let us define now in more detail the GLV framework: Consider a fixed constant 
number of investors (N) (for the extension to a variable number of agents, see 
[18]). At each time step, the wealth of each investor (i) is Wi{t), and the total 
wealth is W{t) = J2iWiit)- 

The time evolution of the Wi{t) is simulated by the following procedure. At 
each time step an investor i is chosen randomly to undergo an event that 
changes its invested wealth Wi{t) — >• Wi{t + 1). The various components in the 
change per unit time dwi/dt are: 

• A deterministic component related to the global status of the economy. 
This term is proportional to the current wealth of the investor through an 
arbitrary coefficient that may depend on time and on all the w's: 

m{wi, ...,WN,t)wi{t) (1) 

One can imagine that the coefficient m{wi, Wn, t) aggregates information 
on economic growth rate, taxes, social benefits, interest rates etc. and it is 
therefore the same for everybody. This is equivalent to the efficient market 
hypothesis [23]. 

• A purely stochastic term which takes into account the specific circumstances 
of each agent. The change of his wealth are still proportional to its currently 
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invested wealth rji{f)wi{f) but the coefficient ry^ is a random number taken 
from a normal distribution with mean < 77 >= and variance < rjf >= Di. 
• There is a social security mechanism, or some fixed relative income that en- 
sures that the investors do not become arbitrarily poor. This term is taken of 
the form Oj bjWj. The coefficients bj represent the relative contribution of 
the individual j to the redistributed wealth (through taxes, donations, pay- 
ments) while flj represent the relative amount that the individual i receives 
from the redistribution (through salaries, services, exchanges, pensions, so- 
cial security). Without loss of generality one may assume J2j bj = 1. In the 
case in which all 6j's are constant: 6j = i/N, the sum w{t) — J2 bjWj reduces 
to the average wealth w{t) — W{t)/N. 

Consequently, in the continuum time limit the GLV dynamics is governed by 
the system of N coupled non-linear differential stochastic equations (in th Ito 
sense) with time dependent coefficients: 

^ = [m{wi..., wn, t) + rii{t)]wi{t) + Oj X] bjWji^) (2) 



Such systems are notoriously difficult to solve or even characterize qualita- 
tively. Yet in the present case, in the limit N ^ oo, (and for positive, not too 
unequal and respectively bi) the probability distribution of relative wealths 
Wi/w is completely under analytic control in spite of the fact that the global 
wealth is very non-stationary and can have arbitrary ups and downs (corre- 
sponding to booms and crashes/ recession). In particular, with the notation, 

« = (3) 
j 

the relative wealth: 

Xi{t) = Wi{t)/w{t) (4) 



has been shown [19-21] to converge even in nonstationary conditions to a 
probability distribution that is m{wi, ...wat, t)-independent: 



Consequently, modulo important finite N corrections [18] which are outside the 
present scope, the dynamics (2) insures in the range Xi > Xmin = 1/(1 + -^^/^^) 
a power law with 

a^l + 2a/D. (6) 
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In the real measurements a has been found to be roughly constant around 1.5 
in the last 100 years in all the western economies [18,20]. 

This value a ~ 3/2 has been related to the intrinsic human biological con- 
straints through the formula a ~ L/(L — 1) where L is the average number of 
dependents on the average wealth [19]. 

In the simulations presented here, we will use a discrete version of Eq. 2 with 
a particular choice of the form of the random factor and of the parameters: 
'm{wi, wn, t) — —a — iJ>w{t), Di — D, Ui — a, hi — 1/N: 

Wiit + 1) = Wi{t) * e'^^W + - Wi{t)) - ^Wi{t)W{t) (7) 

This specific choice oim{w, t) is the minimal form, which embodies the 2 main 
relevant economic facts: 

• the term —awi(t) represents proportional taxation, while 

• the term —-^W{t) models hmiting global factors such as inflation (when 
the total numerary wealth in the system is rising, the real value of each 
investors wealth is decreasing proportionally.). 

The simple Eq. 7 recreates many of the observed features of stock markets, 
mainly the observed distribution of the wealth, the power law in the distri- 
bution of market returns and the long term rise in the total wealth of the 
investors, it fails to reproduce the long term correlations in the volatility of 
the market. We will then show how to amend this problem and how this 
influences the market stabihty. 



3 Volatility correlations in the simple GLV model 

The main focus of the present paper is the market volatility. Let us therefore 
describe first the definition, measurement and properties of this quantity in the 
simple GLV model Eq. 7. We will then study in the next section the modified 
GLV model in which the volatility determines the variance {D) of the random 
factor Tj. 

The market return at time t is defined as: 

The volatility is defined as the average of the square of the returns over a 
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certain time period (we take it as N unit time steps): 

V =<[ln{W{t + l)/W{t))f >N (9) 



The change in W between each time step is small and one can replace eq 9 
by: 



{x,{t) * (e^^W - 1) + a(l - - ^x,{t)W{t)r 

=< x,2>(L> + a^)^^^< xJ>D 



(10) 



One sees that V depends only on the distribution of the relative wealth Xi{t). 

As mentioned above, in GLV the probability distribution of the individual 
relative wealth does not change even in the presence of significant variations in 
the total wealth. Therefore, the classical GLV, the volatility is also a stochastic 
variable with a static distribution. In particular, the volatility inherits the 
scaling properties of the relative wealth distribution. In fact it turns out that 
the volatihty has a Levy distribution with an index of 1 + | (Figure 1). 

The experimentally observed long term correlation of the volatility is a power 
law i.e: < ^^^^^^J'^^ oc t~^, with an exponent close to (5 ~ 0. By contrast, 
the time auto-corrclation of the volatility decays exponentially in the simple 
GLV model. In order to reproduce the experimentally observed property we 
add below an auto-catalytic dependence of D on V(t). 



4 Dynamic volatility 



The volatility represents the "nervousness" of the market. It also measures 
the fraction of money an investor can expect to win or loose by investing in 
the market during a certain time interval. Therefore it is natural to assume 
that the variance D of the random factor 77 (i) in Eq. (7) is in fact a function 
of the volatility V{t): 

D^giV) (11) 
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We will further assume that g{V) can be parameterized by a power of exponent 
n: 

9{V) = ciV^ (12) 

In order to close the feed-back loop one has to estimate the dependence of 
the market volatility V on the variance D of the fluctuations of the individual 

wealth using equation 10. One can measure < > as a function of D (Figure 
2). This dependence turns out to be linear in the range of values used in 
the model (0.03-0.1). Thus in the range of values used in our simulation, one 
obtains that: 

V{D) « C2D'^ (13) 

In fact we measured V{D) for the simple GLV model and verified that it does 
fit this function (Figure 3). 

We are now in the position to estimate the stability of the system as a function 
of n. By using the Eqs 11 and 13 we obtain the iterative equation describing 
the dynamics of the volatility: 

V{t + 1) = C2D{t + 1)' = C2Ciy(i)'" (14) 
The condition for a stationary dynamics ("fix point") is therefore: 

where c — ciC2 

One can now replace Eq 14 with a continuous dynamics: 

V{t + 1) - V{t) = cV(t)(y(t)2"-i - -) 

c 

^T/ = cl/(t)(l^(t)'"-^-^) (16) 
We can now now estimate the stability around the steady state: 

AV = c{Vf, + AV){{Vf, + AVf^-' - -J 

AV = c{2n - 1)V^^-^AV 

AV = (2n - 1)AV (17) 



7 



The fix point Eq. 15 is therefore stable for 2n — 1 < 0, and unstable for 
2n — 1 > 0. If n is exactly n = 1/2 then equation 14 will have a marginal 
steady state at any value of if c = 1. It will diverge for any value of V if 
c > 1, and it will converge if c < 1 to a low enough value of V where the second 
order dynamics will take effect. When the fix point is unstable (n > 1/2) the 
dynamics will lead to zero volatility if we start below the fix point and to 
infinity if we start above the fix point. 

In order to check this effect we simulated the system for various n values: n=l 
(figure 4), n = 1/3 (figure 7,8) and n = 1/2 (figure 5,6) 

• When n = 1 (i.e 2n — 1 > 0) The values of V and D both diverge. The 
divergence mechanism is not driven by the rise in the wealth of the investors, 
or by better investment. It is driven purely by the rise in the variance of the 
market. All the feedback interplay takes place between x (the normalized 
value of the investors wealth) and the volatility V(t). The total W(t) does 
not play an active role in this feedback loop. 

• When n = 1/3. The values of the volatility V and of D both stabilize. This 
will be the case that will be further analyzed. 

• If n = 1/2, and we use c > 1 the volatility diverges, as expected, while 
if c < 1 the values of D and V converge to a very low value. We will 
further investigate the dynamics that leads to this steady state. Interestingly 
enough, this value at the border between divergence and stability seems to 
be favored by the actual market observations. There might exists an self- 
organization argument that explains this fact. 

We assume the diverging case does not represent realistic situations though the 
implied regime might be likened to some of the large fiuctuations experienced 
in the latest years by the Nasdaq index. 

However, we did not consider here, the risk-adversity of the investors might 
be the factor which ultimately prevents the unlimited price increase together 
with the volatility divergence. We also neglected adverse market effects related 
with the largest investors bidding against themselves. 

In the converging cases (n = 1/3 and n — 1/2, c < 1) long-term volatility 
correlations exist (figure 8 and 9). In order to measure the power of the cor- 
relations we made a best-fit estimation to the log of the correlations with the 
log of T and got for this specific set of parameters an exponent of 5 = 0.5. 
Note that this exponent is sensitive to the values of n and c. We plan now to 
study in a future publication the bounds on the parameters n and c that can 
be deduced from the experimental measurments of the volatility corelations. 
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5 Discussion 



The GLV formalism explains elegantly basic features of wealth distribution 
and markets dynamics. One of the main flaws of this formalism was its failure 
to predict the long term corclation in the markets volatility. Long term cor- 
relations imply a long term market memory. Such a memory is not included 
in the GLV formalism, which is based only on the current state of the mar- 
ket. In order for the GLV formalism to include long term corelations memory 
must emerge from the dynamics of the market, instead of being imposed on 
it externally. 

We proposed a simple mechanism that can explain both the long term core- 
lations of the volatility, and the apparent divergence of some markets. We 
propose that the memory is due to a positive feedback between the volatility 
and the nervousness of the traders. The efficient market hypothesis [23] re- 
quires that there can be no direct link between the average gain of the traders 
and any of the measured properties of the market. Thus long term corelations 
cannot be due to a positive feedback on the gains. Thus the next natural 
candidate for this feedback is the standard deviation of the gains. 

The efficiency of the market seems to be also in contradiction with the fast 
rise observed in some markets. We therefore propose here that this rise is not 
due to the average gain of each investor, but to the standard deviation of this 
gain. As a result, one explains the sustained rise in the prices as the result of 
the auto-catalytic feedback of the standard deviation of the traders gain on 
itself. 

We have shown using a stability analysis that a long term rise in the markets 
value may take place for certain values of n the exponent parametrizing the 
dependence of the nervousness of the traders on the volatility. 

If the trade volume increases sharply as a function of the a rise in the volatility 
(very nervous traders) this will lead to a divergent rise in the market prices. 

If on the contrary, the traders are very calm (a weak dependence of the stan- 
dard deviation D of the individual returns on the market volatility V), the 
markets average returns will be stable, and the only trace of this dependence 
will be the long term power-law auto-corelations of the volatility. 

In a previous paper [20] the Constance of the individual relative wealth Pareto 
distribution during the last century was explained in terms of general socio- 
logical factors (size of families...) [19]. In the present paper we showed that 
the long term evolution of the market may be due to psychological factors. 
Even very weak autocatalytic effects of the volatility on the individuals' ex- 
pectations, can determine the long term evolution of the markets. 
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Figure 1 - The volatility distributions. This distribution has a power law tail 
over at least three orders of magnitude. The distribution was measured for a 
constant value of D (the variance of the investors gain) , and after the average 
wealth has stabiliesed due to the non hnear competition factor. The exponent 
of the volatility power distribution is determined by the exponent of the Pareto 
relative wealth distribution. 

Figure 2 - The average of the square of the relative individual wealth < xj >i is 
found (in the relevant parameters range) to depend linearly on the variance in 
the investors random gains D —< rjf >. Deviations from the linear fit occur:- 
at very low values of D, - at large values D ^ 1. For such values of D, most 
of the total wealth is in the hands of the wealthiest agent: x„iax = 1 and all 
other agents have Xi ~ 0. Consequently < > saturates as it approaches its 
maximal possible value 1. However these parameters ranges where linearity is 
violated, are outside the interest of the present paper. 

Figure 3 - The square root of the volatility as a function of the variance in the 
gain (D), and a linear fit. The gain variance was varied and for each variance 
we measured the algebric average of the volatihty once the total wealth reached 
equilibrium. Note that the volatility itself has a very large variance, and these 
results represents the average over a very long time. Thus at small time scales 
the average volatility may be very different than its long term average. 

Figure 4 - The average wealth and volatility, when the variance in the investors 
gain (D) is linearly dependent on the volatility. These represents a very "ner- 
vous" market in which the einvestors are very sensitive to variations in the 
stock values. The dynamics in this case diverge through the positive feedback 
loop between the gain variance and the volatility. 

Figure 5 - The average wealth and volatility, when D is proportional to the 
square root of the volatility. This represents a marginally steady state in which 
the investors sensitivy to the volatility is precisely inverse to the sensitivity 
of the volatihty on the D. The average volatility is stable, however very larfe 
fiuctuations around the average are observed. 

Figure 6 - This simulation is simlar to the one presented in figure 5, with 
a higher value of c. The high value of c leads to divergent total wealth and 
volatility. 

Figure 7 - The simulated volatility with the variance proportional to the cubic 
root of the volatility. In this case the volatihty is stable but it has long term 
corelations. The two drawings represent the volatility fiuctuations at different 
time scales. 

Figure 8 - Long time (r) auto-correlations of the market volatility < V{t + 
T)V{t) > / < V{t)'^ >. The graph was obtained from runs in which the 
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square standard deviation of the individual returns D was proportional to 
V^^^. One observes a straight line fit on the double-logarithmic plot, indicative 
of a power- law auto-correlations decay t~^. 
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Divergence of total wealth 



Time [arbitrary units] 
Divergence of volatiiity 



Time [arbitrary units] 
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Fig. 4. 




Fig. 5. 





Fig. 6. 
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